Scale invariance and viscosity of a two-dimensional Fermi gas 
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We investigate collective excitations of a harmonically trapped two-dimensional Fermi gas from 
the coUisionless (zero sound) to the hydrodynamic (first sound) regime. The breathing mode, which 
is sensitive to the equation of state, is observed with an undamped amplitude at a frequency two 
times the dipole mode frequency for a large range of interaction strengths and different temperatures. 
This provides evidence for a dynamical SO(2, 1) scaling symmetry of the two-dimensional Fermi gas. 
Moreover, we investigate the quadrupole mode to measure the shear viscosity of the two-dimensional 
gas and study its temperature dependence. 
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Scale invariant behaviour plays an important role in 
many branches of physics. It is encountered both in fluc- 
tuations near critical points of phase transitions [1| and 
in particle physics when masses become unimportant at 
high energies j2,]. In non-relativistic quantum mechanics, 
scale invariance implies that the Hamiltonian H{x) scales 
under a dilatation of the spatial coordinate x ^ \x ac- 
cording to H{\x) — ^ H{x)/\^. This scaling symmetry 
allows one to make fundamental statements about ther- 
modynamic properties. For example, even in a strongly 
interacting system in D dimensions, pressure P and en- 
ergy density e are related by the same simple equation of 
state P = 2£/D as an ideal gas. One such system, which 
is of great interest in the ultracold atom, heavy-ion, and 
nuclear astrophysics communities, is the two-component 
Fermi gas in three dimensions interacting via a zero- 
range potential with a scattering length aau. At uni- 
tarity {a^o — ^ oo) it is scale- and conformally-invariant 
with universal properties 0-01 • I^or example, the equa- 
tion of state at zero temperature becomes /i oc Ep, in 
which /i is the chemical potential and Ep is the Fermi 
energy, and the bulk viscosity C vanishes for arbitrary 
temperatures. 

In two-dimensional systems scaling behaviour is more 
subtle. When the two-dimensional scattering length a2D 
^ approaches infinity, the gas becomes non-interacting. 
This implies that at zero temperature one finds /i = 
Ep (X n2D and hence the gas is trivially scale invariant. 
Here, n2D is the density. At finite interaction strength, 
i.e., finite values of a2D, the two-body scattering ampli- 
tude in two dimensions /(g) = — , , 2'^'^ , is momen- 

tum dependent. Evaluating f{q) at a characteristic mo- 
mentum, for example the Fermi wave vector kp, leads to 
a density dependent coupling strength Q . In a quantum 
field theoretical model of the superfiuid Bose gas, it has 
been pointed out that this could give rise to a quantum 
anomaly that breaks scale invariance (loj . 

In the presence of an isotropic harmonic trap, the 
scale invariance of the homogeneous system is replaced 
by a dynamical SO (2,1) scaling symmetry ll[. The 
80(2, 1) group, or "Lorentz" group, is the group of ro- 



tations in 2-1-1 dimensional space-time. For the trapped 
gas, the S0(2, 1) symmetry results in an excitation spec- 
trum with modes spaced by exactly 2w_l 11, [l2|- Here, 
uj±^ denotes the trap frequency of the weakly confined 
axes. This generates a hydrodynamic breathing mode 
at a frequency ujb — independent of the interac- 

tion strength. Moreover, the mode frequency is indepen- 
dent of amplitude and the breathing mode should be un- 
damped. However, the quantum anomaly resulting from 
the density-dependent coupling strength has been pre- 
dicted to shift the hydrodynamic breathing mode of a 
Bose gas by SujbI^b = ^aau/^^ for a^io/lz < 
where denotes the extension of the gas in the strongly 
confined direction. 

In this manuscript, we study an interacting two- 
dimensional Fermi gas using collective modes to inves- 
tigate scale invariance and viscosity. We tune the ratio 
o-3d/Iz between —3 and 0, which provides us with access 
to the hydrodynamic and the coUisionless regimes. Pre- 
vious experiments using ultracold atomic gases to study 
scale invariance of a two-dimensional Hamiltonian were 
limited to weakly interacting bosons in a regime where 



< a^D /Iz ^ 1 U3, UM s-iid to weakly interacting three- 
dimensional gases with a highly elongated symmetry jl6} . 

In our experiment 17,[l8|, we create two-dimensional 
Fermi gases of ^"K atoms in a 50/50 mixture of the 
\F = 9/2, mp = -9/2) and |F = 9/2, -7/2) 
states of the hyperfine ground state manifold. The Fermi 
gas is loaded into the standing wave potential of an op- 
tical lattice to create an array of two-dimensional gases. 
The trapping frequency along the strongly confined di- 
rection is cj^ = 27r X 78 kHz and the radial trapping fre- 
quency is = ^ijjx'^y — 27r X 125 Hz with an anisotropy 
e = \ujx — ijJy\/'2,ijj±_ below 2%. Along the axial direction 
we populate approximately 30 layers of the optical lattice 
potential with an inhomogeneous peak density distribu- 
tion of typically 2 x 10'^ atoms per spin state per 2D gas at 
the center. We tune the interactions by applying a mag- 
netic field close to the Feshbach resonance at 202.15 G 

First, we study the quadrupole mode of the two- 
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FIG. 1. (Color online) a: Frequency 
of the quadrupole mode at T/Tf ~ 
0.47 with Ef ^ hx (8.2 ± 0.7) kHz. 
The dotted line shows the theoret- 
ical prediction of [l^] in the coUi- 
sionless limit and the dashed line at 
cjq/ojx = V2 is the hydrodynamic 
case, b: Damping of the quadrupole 
mode. The solid line shows the fit for 
To (zero sound) and the dashed line 
for Fi (first sound) - see text. The 
dash-dotted line is the damping rate 
of the dipole mode, c: Frequency of 
the breathing mode. For the strong 
excitation we have T/Tp ~ 0.37 and 
EF = hx (5.4 ± 0.8) kHz and for the 
weak excitation T/Tf = 0.42 and 
Ef ^ hx (8.1 ± 1.1) kHz. d: Damp- 
ing of the breathing mode. The 
dash-dotted line is the damping rate 
of the dipole mode. 



dimensional Fermi gas in order to identify the coUisionless 
and the hydrodynamic regimes. The quadrupole mode 
has the incompressible velocity field VQ(r) = b[xex — 
yBy] cos(wgi) with a constant b, and corresponds to a 
surface mode of the gas. In the coUisionless limit, sur- 
face modes are analogous to zero sound. At zero tem- 
perature the quadrupole mode frequency is predicted to 
be u}Q = \/2{2 - g)/{l - g)uj± with g = l/ln(fci?a2D) 
[l3| . A collective hydrodynamic mode, in contrast, cor- 
responds to a first sound mode. The frequency of the 
hydrodynamic quadrupole mode ujq — y/2uj± 20, 21 1 is 
independent of the equation of state because its incom- 
pressible flow pattern prevents a change of the internal 
energy during the oscillation [2^ . 

We excite the quadrupole mode by adiabatically intro- 
ducing a small anisotropy to the two-dimensional har- 
monic oscillator potential using additional laser beams 
while maintaining x « const, and then abruptly 
returning to the original trapping configuration. The 
atomic cloud oscillates freely in this potential for up to 
14 ms until we switch it off and take an absorption image 
after 12 ms of time of flight. The velocity amplitude of 
the excitation is 10% of the Fermi velocity vp — hkp/m. 
We determine the radius of the cloud in the x- and the 
y-direction and fit their difference to measure loq p3| . 
Owing to the change in gravitational sag during the ex- 
citation, we simultaneously excite small dipole (center- 
of-mass) oscillations primarily in the vertical y-direction. 
We use these dipole oscillations to calibrate Lo^^y, of which 
ijjy has the smaller error because of the larger oscilla- 
tion amplitude. The decay rate of the bare dipole mode 
To/^i_ = 0.04 ± 0.01 is most probably caused by the 
weak anharmonicity of our Gaussian trapping potential. 

In Figure [l}i we show the quadrupole mode frequency 
UJQ of the two-dimensional Fermi gas. For large values 
of the interaction parameter ln(fci?a2D) we are in the 



coUisionless regime and observe loq ~ 2aj_L, in agree- 
ment with the theoretical expectation [l3| . As we in- 
crease the interaction strength, i.e., decrease the value 
of ln(fci?a2_D), we enter the hydrodynamic regime. This 
transition is marked by a sharp decrease of the collective 
mode frequency to ojq ~ \/2ajj^. Theoretically, we ex- 
pect the transition from the coUisionless to the hydrody- 
namic regime when the collision rate 70 equals the mode 
frequency ujq. We estimate 70 = — /in2Dlm[/(fci?)]/m 
using the optical theorem for the scattering cross sec- 
tion a = — Im[/((7)]/g. For very small deviations from 
the equilibrium distribution, the collision rate 70 is sup- 
pressed by a factor (T jTp)"^ < 1 owing to the restriction 
of phase space because of Pauli's exclusion principle. Us- 
ing our average density n2D ~ 6 x lO^^m"^ and temper- 
ature T/Tp — 0.47 we estimate the transition from the 
hydrodynamic to the coUisionless regime at an interac- 
tion parameter lii{kpa2D) ~ 3, in good agreement with 
the observed mode frequency change. 

In Figure [T]d we show the damping rate Tq of the 
quadrupole mode. The zero sound mode in the coUi- 
sionless regime is damped by collisions which disrupt the 
coherent quasiparticle motion. Hence, the mode damp- 
ing rate Fq scales proportional to the normalized col- 
lision rate ^q/ujq with an asymptotic behaviour Fq oc 
[\n{kpa2D)] ^, shown as the solid line. The damping 
rate reaches a constant offset of 0.1 uj±, even for the non- 
interacting gas, which we will investigate in more detail 
below. For the first sound mode in the hydrodynamic 
regime the situation is opposite: collisions are necessary 
to establish local equilibrium and first sound is damped 
by the deviation from this, hence Fi oc ujq/jq, which 
asymptotically is Fi cx l+[21n(/i;i?a2_D)/7r]^ (dashed line). 
In the fit of the proportionality constant we have used 
the same constant offset as for the non-interacting gas. 
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In between the two extremes, the damping rate peaks at 
the transition from the colhsionless to the hydrodynamic 
regime. 

Having identified the hydrodynamic and colhsionless 
regimes, we now turn our attention to the breathing 
mode and the question of scale invariance of the two- 
dimensional Fermi gas. The velocity field of the breath- 
ing mode is vs(r) = b[xex + yCy] cos(w_Bt) with a con- 
stant b. The breathing mode is excited by adiabatically 
decreasing the strength of the two-dimensional confine- 
ment and then abruptly returning to the original trap- 
ping configuration. After an oscillation time of up to 
20 ms we switch off the confinement and take an absorp- 
tion image after 12ms time of flight. We identify the 
breathing mode frequency (Figure [TJ;) by studying the 
cloud radii in x- and y-direction as a function of hold 
time and define ujb = ,Ji^B,x^B,y Again, we use the 
simultaneously excited dipole mode for continuous ref- 
erencing. We study two sets of data of the breathing 
mode as a function of interaction str eng th which differ 
by a factor of 2 in excitation strength [2J] . The weak ex- 
citation strength corresponds to 12% modulation of the 
width after time of fiight. We observe that the mode 
frequency is approximately constant for all interaction 
strengths and averages at ws/wy = 2.00 ± 0.03 for the 
weak excitation and at ws/wy = 1.96±0.01 for the strong 
excitation. We also have not observed any change of the 
mode frequency with the temperature of the gas in the 
range 0.37 < T/Tp < 0.9. 

The independence of the breathing mode frequency 
from interaction strength, oscillation amplitude, and 
temperature and the fact that the mode frequency is very 
close to ujb = 2uj±^ suggest that the gas exhibits a dynam- 
ical S0(2, 1) scaling symmetry. In the hydrodynamic nor- 
mal regime, our result implies that the equation of state is 
polytropic and fi oc n2D ■ Interestingly, this result is com- 
patible with BCS mean field theory at finite temperature 
[2ll [25I . A more precise determination of the equation 
of state, which also takes into account logarithmic inter- 
action energy shifts and other beyond mean-field effects, 
could be obtained from Quantum-Monte Carlo calcula- 
tions, which so far are only available for superfluids at 
zero temperature [2^ . 

The damping of the breathing mode (Figure[T}l) differs 
fundamentally from the quadrupole mode, because it is 
very small and independent of the interaction strength. 
In order to understand the undamped behaviour of the 
breathing mode, we employ a simplified model based on 
classical hydrodynamics. Energy dissipation of the ve- 
locity field VB(r), and hence damping of the coherent 
motion of particles, is caused by the viscosity of the 
gas. The energy dissipation rate E follows from the two- 
dimensional stress tensor [2^ E ~ —\ J (Pr ri(r){dkVi + 
diVk — 6ikV ■ v)^ — / d'^r C(r)(V • v)^ with the shear vis- 
cosity ry(r) and the bulk viscosity C(i")- The time av- 
eraged energy dissipation rate is {E)t = —26^/ (Pr((r), 
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FIG. 2. Difference between the breathing mode frequencies 
of the two different axes at the collisionless-hydrodynamic 
crossover in a slightly anisotropic trap. The insets show the 
correlation plot of the cloud widths along the two axes. 



entirely determined by the bulk viscosity. The amplitude 
damping of the mode is F = {E)t/2{E)t, using the time- 
averaged mechanical energy {E)t = / (Pr r^n{r). In 
our data, the damping rate of the breathing mode av- 
erages near Ts/ujy ~ 0.05, equal to the damping of the 
dipole mode which is dominated by technical limitations 
(e.g. dephasing due to trap anharmonicities) rather than 
by viscous forces. Additionally, we observe the absence 
of damping for different excitation amplitudes and tem- 
peratures. 

As a side remark, we observe the transition from the 
colhsionless to the hydrodynamic regime also directly 
for the breathing mode. In a slightly anisotropic trap 
(e < 0.02) in the colhsionless regime, the oscillation fre- 
quencies of the widths of the cloud are split by 2ewB, 
corresponding to two independent modes (see Figure [2]) . 
When the interaction strength is increased into the hy- 
drodynamic regime the two mode frequencies lock to- 
gether and the whole cloud undergoes a collective breath- 
ing mode at a single frequency. The occurrence of this 
"mode-locking" coincides with the values of ln(fci?a2_D), 
where we observe the collisionless-hydrodynamic transi- 
tion of the quadrupole mode. As insets we show cor- 
relation plots between the widths of the cloud demon- 
strating the change from nearly uncorrelated into highly 
correlated motion. 

We now turn our attention to the temperature depen- 
dent damping of the quadrupole mode. The viscosity of 
a strongly interacting gas plays an important role in the 
investigation of the ratio of viscosity to entropy density 
predicted by the anti-de Sitter/Conformal Field Theory 
correspondence, which so far has been primarily inves- 
tigated in three dimensions p8l - l30| . In two dimensions, 
the shear viscosity of the Fermi gas and its temperature 
dependence have been theoretically investigated only in 
specific limits 



311 . |32| , none of which correspond to our 



experimental situation. Our data for the temperature de- 



4 



1.0 



E 

(0 

Q 



0.0 



S 2.7 
O 5.3 
O 9.7 
O 18 
O 42 
545 



0.0 0.2 0.4 0.6 0.8 1.0 
T/T, 



24 



S 
3 




£ 0.4 

c 

o 

CL 
X 

LU 0.0 



10 250 500 

ln(kpa^J 



10 250 500 



FIG. 3. (Color online) Temperature dependent damping of 
the quadrupole mode, a: Damping rate as a function of 
temperature for various interaction strengths. Ep/h for the 
data sets are (from lowest to highest temperature): 6.4kHz, 
8.2 kHz, 9.0 kHz, and 9.1kHz. Inset: Derived values of 
aiT/Tp) and power-law fits, b: Viscosity amplitude ao- c: 
Power law exponent (3. 



pendence of the damping rate Fg are displayed in Figure 
for various interaction strengths. We use the same hy- 
drodynamical model as above in order to link the shear 
viscosity to the damping rate, which neglects possible 
temperature gradients as well as quantum statistics. The 
time averaged dissipation rate of the quadrupole mode 
is {E)t = —2b^ J cPrr]{r). Using the parametrization 
ri{r) = hn2D{^)o!-{T /Tp) with a dimensionless viscosity 

m{r'^)VQ/2h. We 



aiT/Tp) I33|, we obtain aiT/Tp) 

determine the rms-radius of the cloud -y/ (r^) numerically 
for the non-interacting gas. We fit our data with a power 
law a{T/Tp) = ao x (T/Tp)^ , which is known to pro- 
vide the high-temperature scaling in three dimensions. 
The extracted amplitudes ao and exponents /3 are shown 
in Figures [3Jd and [St, respectively. For \n{kpa2D) > 10, 
i.e. deep in the collisionless regime, we observe no signif- 
icant temperature dependence (i.e. /3 ~ 0) and hence a 
constant a, which at least in part could be due to tech- 
nical limitations such as anharmonicities. In the weakly 
interacting regime, the temperature dependence of the 
viscosity is significant with /3 ~ 1/2. In contrast to three 
dimensions, where the viscosity has been investigated in- 
tensively 0, [28l - l30j , no theoretical prediction exists yet 
for this parameter range of a two-dimensional Fermi gas. 



In conclusion, we have studied collective oscillations of 
a two-dimensional trapped Fermi gas in the collisionless 
and the hydrodynamic regimes. We have observed the 
existence of a breathing mode at two times the trap fre- 
quency, which is invariant against interaction strength, 
amplitude of the excitation, and temperature. More- 
over, this breathing mode is undamped as compared 
to the dipole mode. These observations suggest a dy- 
namic S0(2, 1) scaling symmetry of the trapped two- 
dimensional Fermi gas. In our parameter range we do 
not observe indications for a quantum anomaly breaking 
scale invariance. Using the quadrupole mode, we have 
additionally studied the temperature-dependence of the 
shear viscosity of the two-dimensional Fermi gas. 
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